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1. Introduction 

The helium-hydride ion HeH + is the simplest heteronuclear two-electron system and the 
only well bound molecular system of hydrogen and helium, the most abundant elements 
in the universe. This motivated many calculations and experiments performed on this 
molecular ion which is of interest for the chemistry of astrophysical objects, for the 
tritium neutrino mass experiments, or simply by itself as a model system (see, e.g., 
PiElEllIlElElITJIElEl ED] and references therein). More recently, HeH + is drawing 
special attention due to an experiment performed at the free-electron laser FLASH [IT] 
where the absolute photodissociation cross section at a photon energy of 38.7 eV was 
measured. 

The photoabsorption spectrum of HeH + was previously investigated theoretically, 
together with photoionization, only for the parallel orientation of the molecular axis 
with respect to the field [Sj. A more recent study on the photoionization of HeH + [12] 
treated also the perpendicular orientation, but for the photodissociation into electronic 
bound states, to the authors' knowledge, the only work considering both the parallel and 
the perpendicular orientation is still the one of Basu and Barua from 1984 [7]. There, 
however, the IT contribution taken into account is only the one of the lowest state, while 
the FLASH experiment shows an important contribution of higher lying n states [TT] . 

Motivated by the FLASH experiment, the photodissociation cross section of HeH + 
has been calculated for parallel, perpendicular, as well as isotropic orientations of the 
molecular axis with respect to the field, for transitions starting from the lowest lying 
vibrational state and the first or second rotational level, as well as for a mixture of initial 
vibrational states. Since the experiment measures only the dissociation into He + H + , a 
comparison to the experimental data required resolving different dissociation channels. 
The influence of non-adiabatic effects on the channel-resolved results is discussed. The 
paper is organized as follows: section [2] presents the method and the computational 
details, section [3] discusses the results and a summary is given in Section HI 

2. Method and computational details 

In atomic units (I Eh = 27.21139 eV and 1 ao = 5.29177 • 10~ n m) and the non- 
relativistic dipole approximation the photodissociation cross section is given according 
to Fermi's Golden Rule by 

475-2 

0"M = 9 & K^n'e'J'm', \^>\^ nujm ,)? ' S{uJ ~ [E n , e ,j, m i - E wJmj \) (1) 

C 

where uo is the photon energy and D is the dipole operator in length gauge. For a 
diatomic molecule, it is sufficient to calculate the single-photon absorption cross section 
for a parallel and a perpendicular orientation of the molecular axis with respect to 
the electric field component. The partial parallel (perpendicular) cross section for an 
isotropic molecular ensemble is then obtained by a multiplication of the cross section 
for a molecule oriented parallel (perpendicular) to the field with the statistical factor 
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g = 1/3 (g = 2/3). The total cross section for an isotropic molecular ensemble is simply 
the sum of these two partial cross sections (see also [H [13] for more details on g). 

In equation (JTJ, ^ nu jmj and ^n'e'j'm'j are the rovibronic wavefunctions of the initial 
and final state, with energies E nvJmj and E n i e rji m / respectively. Within the Born- 
Oppenheimer (BO) approximation the former is given as 

^n V JmAr Xl r 2l R) = Mrur 2 ;R) -^-L Y Jmj (p) (2) 

where ip n {fiiT 2 ;R) is the electronic, $>„](R)/R the vibrational, and the spherical 
harmonic Yj mj the rotational wavefunction. The index n specifies the electronic state 
(BO potential curve) including its symmetry (the only relevant ones in this work are 
X S and 1 I1), while v and {J,mj} are the vibrational and rotational quantum numbers, 
respectively. The final-state wavefunction can be written in an analogous way, but in 
the here studied case of photodissociation e' is a continuous index, while n, n', u, J, and 
J' are discrete ones. Correspondingly, ip n , ip n > } , and, of course, Yj mj as well as K/' m ^ 
are normalized to unity, while &™)> is energy normalized, = 5(e' — e"). 

The transition amplitudes in equation ([TJ) are obtained by a calculation of the 
integral 

(Vn'e'J'm'plVnvJmj) = J (R) D n , n (R) (R) dR (3) 

where D n i n (R) denotes the electronic transition dipole moment as a function of the 
internuclear distance. It is defined by 

D n , n {R) = J r 2 ; R) d ^„(n, f 2 ; R) d 3 r (4) 

where d is the electronic dipole operator in length gauge, d = e ■ (r\ + r 2 ) with e being 
the polarization vector. (The contribution of the nuclear dipole operator vanishes, since 
in the present work only electronic transitions with n ^ n' are considered.) 

The electronic wavefunctions and transition dipole moments are calculated using 
the method described in [2] which is thus only briefly sketched in the following. First, 
the one-electron Schrodinger equation (OESE) is solved in an elliptical "box" within 
the fixed-nuclei approximation, using a 5-spline basis set and the prolate-spheroidal 
coordinates (1 < £ < oo, — 1 < rj < 1, < < 2n). The box boundary is defined 
by £max- The two-electron Schrodinger equation is solved using the configuration- 
interaction (CI) method in which the configurations are expanded in terms of the one- 
electron orbitals obtained from the solution of the OESE. This yields the BO potential 
curves and the two-electron wavefunctions ip. The latter are used to evaluate the 
electronic transition dipole moments D n / n . All the results shown in this work were 
calculated in the length gauge, but agreement with velocity gauge was also investigated. 

The nuclear motion is solved in the adiabatic BO potential curves by expanding 
the nuclear wavefunctions in B splines for the radial part times spherical harmonics 
describing the angular part. The radial Schrodinger equation is solved in a spherical 
"box" whose boundary is defined by R max . The box leads to a discretization of the 
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vibrational continua. The resulting vibrational wavefunctions are therefore all 

discrete and normalized to unity. For the dissociative continuum the renormalization 
factors that lead to an energy normalization can simply be obtained from the density 
of states. 

The present calculations were directly motivated by the recent FEL experiment 
Therefore, the basis sets were chosen for the purpose of computing photodissociation 
cross sections that take into account the contribution of as many Rydberg states as 
possible (the experimental results indicate the importance of highly excited states). In 
the present calculation 35 X and 25 II electronic bound states were finally included. A 
physical box (r max = -|£ max ) equal to 40 a was used, changing the size of the elliptical 
box £ max synchronously with the variation of the internuclear distance R, 40 B splines 
of order 8 were used for the £ direction and 20 B splines of order 6 for the r] direction. 
The BO potential curves were calculated in an R interval in between 0.1 and 60 a . The 
latter limit was mainly needed for a determination of the dissociation channel, i. e. to 
find out into which atomic or ionic fragments HeH + dissociates. The basis set used for 
calculating the vibrational wavefunctions $ comprised a box R max = 20 a and 900 B 
splines of order 10. 

3. Results 

Calculations were performed for transitions starting from the lowest lying vibrational 
state of the electronic ground state to all the calculated electronic bound states, for 
both parallel and perpendicular orientation of the molecular axis with respect to the 
field. The results shown here are for transitions starting in the rotational level J = 1, 
which is the most probable value of J at room temperature [7], and going to J' = 0, 2 
for the parallel orientation (according to the selection rule AJ = ±1, AJ ^ 0) and to 
J' = 0, 1, 2 for the perpendicular transitions (selection rule: AJ = 0, ± 1). Since previous 
calculations [8] were performed for transitions starting from J = 0, the spectrum for 
these transitions is also shown for comparison, but the final result is not significantly 
influenced by the rotational motion within the considered J values. 

3.1. Energies, cross section, and comparison to previous results 

Calculations of electronic energies and dipole moments for different values of the 
internuclear separation were performed. For checking the accuracy of these calculations, 
results were compared to literature (|2l (3j [8]) and good agreement was found for the 
basis sets described in Section [2j Plots of the computed potential curves for the X S 
and 1 I1 symmetries are shown in figures [1] and [21 respectively. These potential curves 
were used in the subsequent calculation of the photodissociation cross sections. Clearly, 
all potential curves are purely dissociative within the Franck-Condon window of the 
electronic ground state. One-photon transitions from the ground to the excited states 
lead therefore to dissociation. 
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Figure 1. The potential curves of the 1 S states of HeH+ (included in the present 
calculation): the full black lines represent states dissociating adiabatically into He + 
H+, while the red dashes represent the H + He"'" channel. The black arrow points at 
an (also encircled) avoided crossing that is expected to influence the channel-resolved 
results. The blue vertical arrow marks the experimental photon energy used in 

The gauge dependence of the results was investigated by comparing the electronic 
transition dipole moments obtained in either the length or the velocity gauge. The 
latter moments are obtained from equation (TjJ, if the dipole operator is expressed as 
d v = e ■ (Vi + V2) and multiplied by - (u> being the photon energy) for comparing to 
the length gauge results. Figure [3] shows a comparison of the transition dipole moments 
obtained in the two gauges for the equilibrium internuclear separation R eq = 1.45 ao and 
for a parallel orientation (transitions to states). The overall agreement is evidently 
very good. Noticeable deviations are only found for some very high lying electronic 
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Figure 2. As figure Q] but the 1 I1 states are shown. 



Rydberg states that possess, however, a very small transition moment. These deviations 
do therefore have an extremely small influence on the cross-sections obtained from all 
electronic states. 

In figure H] the calculated photodissociation cross section for a parallel alignment of 
the molecular axis with respect to the photon polarization vector (g = 1) is compared 
to a previous result [8] in which explicitly correlated basis functions had been used 
for the evaluation of the electronic wavef unctions. Both results were obtained for 
transitions starting from the absolute rovibronic ground state of the HeH + molecule 
[y = 0, J = 0). The most evident difference between the two spectra lies in the 
position and height of their maxima. This difference is due to the fact that the electronic 
wavefunctions especially for the ground and the first excited states are more accurately 
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Figure 3. The absolute value of electronic transition dipole moments for the 
equilibrium internuclear distance R = 1.45 ao, parallel orientation. The black circles 
are the length gauge results, while the red crosses are the velocity gauge dipole 
moments. 



represented in [8] than in the present calculation, as is concluded from a comparison of 
the corresponding energies. 

In references |14[ [15] it was demonstrated that the present approach can yield 
extremely accurate two-electron wavef unctions, but this requires a judicious choice of the 
basis set, if the size of the calculation should be kept reasonably small. Most importantly, 
different electronic states possess different optimal basis sets. On the other hand, the 
goal of the present work was to consider as many excited bound states as possible and 
to calculate all of them with sufficient and comparable accuracy. The most convenient 
way to do this is to use a single basis set for all states, since it avoids non-orthogonality 
problems and preserves sum-rule completeness. This leads to a compromise for the 
results for the ground and first excited states. It was checked that by improving the 
ground and first excited states the agreement to [8] in the peak region is also improving. 
The remaining part of the spectrum that stems from transitions to higher-lying electronic 
states is, however, in good agreement with the previous result in [8]. As is discussed 
below, for a comparison to the experiment reported in [11] the accuracy of the first 
excited electronic state plays practically no role and thus the good agreement with the 
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Figure 4. Comparison of the results obtained for the parallel orientation of the 
molecular axis with respect to the field to the results from [8]: the full black line is 
the result of the present calculation for J = 1 — > J' = 0, 2, the blue chain is present 
calculation for J = — > J' = 1, and the red dashed line is the result from [8] (obtained 
foJ = fl-+J' = l). 

previous calculation for the higher photon energies is more important for the purpose 
of this work. 

Figure H] shows also the parallel photodissociation cross section for the first excited 
rotational state (J = 1) which, according to |7J, is the most probable rotational level 
at room temperature. Accidentally, the result for J = 1 agrees better to the previous 
theoretical calculation performed for J = 0. From this accidental agreement and the 
discussion above it is clear that the influence of rotational excitation is in general rather 
small and most prominent for the transition to the first excited electronic state. It was, 
however, checked that the results shown in the following are all relatively insensitive to 
the initial rotational state, as long as J is not very large. 

Figure [5] shows the photodissociation cross section of HeH + for an isotropic 
orientation of the molecular axis with respect to the field, as well as the parallel (g — 1/3) 
and perpendicular (g = 2/3) contributions to it, for transitions starting from the lowest 
lying vibrational state v = and the first excited rotational level J = 1. It can be 
seen in figure [5] that until about 27.5 eV the spectrum is practically identical to the 
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Figure 5. Total photodissociation cross section of HeH+ after averaging over 
an isotropic orientational distribution (blue solid). Also shown are the partial 
contributions from parallel (black dashes) and perpendicular (red chain) orientations. 
The calculation is performed for transitions starting from the lowest lying vibrational 
state v — and from J = l.The crosses are the theoretical result from [7] and the 
vertical chain at 38.7 eV marks the experimental photon energy. 

parallel contribution, while above about 29 eV the perpendicular contribution is clearly 
dominant. Similarly to the parallel case the perpendicular spectrum shows a dominant 
low-energy peak followed by a rather structureless tail. In fact, the tail looks more like 
a second peak, but it is composed of a number of electronic transitions. In comparison 
to the tail for a parallel alignment the one in the perpendicular case is confined to a 
smaller energy window. 

An older calculation [7J (performed using the length gauge) took into account only 
three excited states of HeH + , namely the first two excited E states and the first IT 
state. Figure [5] shows that for these particular states the agreement between the present 
calculation and the one from [7J is excellent, but at the photon energy used in the 
FLASH experiment higher lying states clearly dominate the spectrum. 
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3.2. Comparison to experiment: channel separation, diabatic effects, and vibrational 
excitation 

The heteronuclear HeH+ ion can dissociate into two different channels: He + H + and 
H + He + . The FLASH experiment [11] measured only the cross section given by the 
dissociation of HeH + into He + H + . A comparison requires thus a separation of the two 
channels. Within the adiabatic approximation the asymptotic long-range behaviour, 
i. e. the separated atom limit at R — ■> oo, defines the channel into which the population 
of a given adiabatic electronic state dissociates. At R = 60 a the potential curves were 
found to have almost reached the separated atom limit. A comparison of the asymptotic 
potential-curve energies with the sum of the energies of an independent atom and ion 
(He + H + or H + He + ) in their various states yielded a well-defined channel assignment. 
The result is shown in figures [1] and El Most of the lower lying states dissociate into 
He + H+ while the highest lying Rydberg states dissociate into H + He + . The two 
blocks are rather well separated, except for three states (for both X S and 1 U symmetry) 
that dissociate into H + He + but lie in the energy regime dominated by the He + H + 
channel. 

After considering only transitions to electronic states that are assigned to the He 
+ H + channel the partial photodissociation cross section shown in figure [6] is obtained. 
The result (obtained for v = and J = 1) shows in comparison to figure [5] that overall 
dissociation into the He + H + channel is clearly less favoured. However, the relative 
importance of the two channels varies strongly with photon energy. Below about 27 eV 
the spectrum stems solely from the first excited state of X S symmetry which dissociates 
into H + He+. In between about 30 and 35 eV the total dissociation cross section is 
dominated by the contribution from the lowest lying x n state that also dissociates into 
H + He + . A smaller contribution comes, however, from the second excited X E state. 
This state is responsible for the first peak observed in the He + H + channel in figure 
|6j Thus the parallel contribution is equal to the total cross section of this channel in 
that energy range. Above about 35 eV many other states contribute, and the cross 
sections for the two dissociation channels are almost of equal size, as a comparison of 
figures [5] and [6] shows. Since the majority of the states dissociating into He + H + lies 
energetically above the ones dissociating into H + He + , the maximum of the second 
peak of the He + H+ cross section is shifted to slightly higher energy than the one of 
the total dissociation cross section. 

For comparison, figure [6] shows also the experimentally determined absolute partial 
dissociation cross section at 38.7 eV for the He + H+ channel as well as its parallel 
contribution [TT]. The agreement is almost perfect, but further aspects have to be 
considered for a direct comparison. First, the present analysis is based on an adiabatic 
model. Second, it appears to be clear that the HeH + ions investigated in the experiment 
were not in their vibrational ground state. These two issues are thus discussed in the 
remainder of this section. 

The potential curves in figures [Hand [2] show a number of avoided crossings. Clearly, 
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Figure 6. Partial photodissociation cross section of HcH+ yielding He + H + . The 
calculations are performed for transitions starting from the lowest lying vibrational 
state v = and from J = 1 . The blue square is the experimental total cross section and 
the black circle is the experimental parallel contribution (both from [11]). Adiabatic 
limit: the blue solid line represents the isotropic orientation of the molecular axis 
with respect to the field. The parallel (black chain) and perpendicular (red dots) 
contributions are also shown. Diabatic limit: the magenta dashes represent the 
isotropic orientation, while the green chain is the parallel contribution. 



at these crossings there is some probability for a transition from one adiabatic curve 
to the other. In the present context not all avoided crossings are, however, of equal 
importance. Those between states dissociating into the same dissociation channel do 
not change the partial cross section into, e.g., the He + H + channel. (Clearly, these 
avoided crossing have to be considered in a calculation of other differential spectra 
like the kinetic-energy distribution of the ionic or atomic fragments.) A more careful 
analysis that takes also the relative importance of the different electronic states (at the 
experimental photon energy) into account shows that the most critical avoided crossing 
for the present discussion is the one occuring close to R = 4.0 between the 2nd and 
3rd excited state of X S symmetry (see figured]). (This avoided crossing is already visible 
in the theoretical data in [2] [16] and thus reproduced by other numerical methods.) 
A standard procedure for estimating the importance of diabatic effects at an avoided 



Photodissociation of the HeH + molecular ion 



12 



curve crossing is provided by the Landau-Zener approximation (see, e.g., Chapter 10 
in [T7] and references therein). Within this model the probability for a transition from 
one adiabatic curve to the other is given by 



where R c is the internuclear separation at which the avoided crossing occurs, Hi2(Rc) is 
the diabatic coupling which can be calculated from the relation E 2 — E 1 = 2H 12 (Ei and 
E 2 being the adiabatic energies), the velocity v is obtained from the kinetic energy and 
F 12 are the slopes of the diabatic curves at R c . In the present case, it turns out that 
the Landau-Zener formula is difficult to apply, since the two curves cross in a region 
where their behaviour is not at all linear. (A removal of the nuclear repulsion term does 
not resolve this problem.) As a consequence, it is difficult to estimate the slopes and 
the result of the Landau-Zener estimate has such a large error bar that it is practically 
needless. Consequently, figure M shows the cross sections in the two extreme cases: the 
already discussed adiabatic limit (P" 2 = 0) and the fully diabatic one (P" 2 = 1). From 
the figure one may conclude that the influence of possible diabatic effects on the cross 
section at the photon energy used in the experiment is not very relevant. The most 
noticeable changes would take place at lower photon energies. 

All the spectra shown so far were obtained by assuming that the system starts in 
the lowest lying vibrational state v — 0. This is valid for thermal distributions at not too 
high temperatures. However, due to the way the HeH + ions are generated in the FLASH 
experiment reported in [11] the detailed vibrational distribution is unknown, but not 
expected to be thermal. In a more recent experiment on He^ in which the ions were 
produced in the same way (by heavy particle collisions) the vibrational distribution was 
determined [18] . In this case the initial-state distribution was found to be distributed 
over at least 5 vibrational states iy = to 4). Since the binding energy of He^ is 
similar to that of HeH + , one may expect a similar vibrational-state distribution for 
both ions. Therefore, photodissociation cross sections for transitions starting from 
different vibrational states were computed. The ones corresponding to the initial states 
u = 0,1, 2, and 3 are shown in figure El 

At the experimental photon energy the results for both the adiabatic and diabatic 
limit are very close to each other, but there are evident differences at lower photon 
energies. Accidentally, the photon energy in [TT] probes the spectrum of the HeH + ion 
at a point where the cross sections of the v = and v = 1 states cross each other. 
As a consequence, the performed experiment is rather insensitive to a redistribution 
of population between the states with v = and v = 1. The population of higher 
vibrational states iy > 1) leads on the other hand to a decrease of the cross section at 
38.7 eV. The trend visible for v = 2 and v = 3 continues for higher values of v, as was 
checked though not explicitly shown. 




(5) 
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Figure 7. Partial photodissociation cross section of HeH+ decaying into He + H + 
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Figure 8. Photodissociation of HcH+ into He + H + - comparison between theoretical 
calculations (adiabatic - red dashes, and diabatic estimate - black solid line) for 
transitions starting from a mixture of vibrational levels (33% v = 0, 45% v = 1, 
and 22% v > 2) and initial rotational level J = 1 with the experimental (black square) 
result UTI. 



If one assumes the initial vibrational distribution of HeH + in [TT] to be identical to 
the one of He^ 4 " in [18] . the photodissociation cross section shown in figure M is obtained. 
The vibrational distribution is in this case modelled in the following way: 33% of the 
vibrational population was considered to be in the v = level, 45% in v — 1, and 
22 % are equally distributed over the states v = 2,3,4. Compared to figure [6] the 
influence of the vibrational distribution is obvious: the shape of the spectrum changes 
and the value of the cross section at the experimental photon energy (38.7 eV) increases 
slightly (from a ~ 1.50 Mb - see figure [B]- to a ~ 1.70 Mb). Thus the agreement to 
experiment decreases, but remains still well within the error bars of the measured cross 
section. In view of the lack of knowledge about the experimental initial vibrational 
distribution, it is, however, impossible to finally conclude on the degree of agreement 
between theory and experiment. As the present study shows, this distribution has a 
rather strong influence on the photodissociation cross section. In fact, the theoretical 
results may indicate a narrower vibrational distribution than the one found for He2 + in 
[T8] that was assumed to be similar to the one of HeH + in [11] . It would certainly be of 
interest to have experimental data also for other photon energies. This may also help 
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to (indirectly) determine the initial vibrational-state distribution. Furthermore, these 
data may be an interesting direct measure of the adiabaticity at the discussed avoided 
crossing between the 2nd and 3rd excited state of X S symmetry. 

4. Summary 

In this work the photodissociation cross section of HcH+ was computed. A previous 
theoretical result that considered a parallel alignment of the molecular axis was 
confirmed. A calculation for a perpendicular alignment was performed that considered 
a large number of electronic states in order to obtain a converged photodissociation 
cross section. In order to compare to a recent experiment, the partial cross sections for 
dissociation into either He + H + or H + He + were determined. This was done within 
the adiabatic approximation, but possible diabatic effects were estimated. Different 
initial vibrational and rotational levels were considered and photodissociation spectra 
for transitions starting from the rotational levels J = 0, 1 and vibrational levels 
v = 0,1,2,3,4 were presented. The influence of the vibrational motion on the final 
spectra was investigated and compared to the experimental findings. Agreement within 
the experimental error bar was found, but it was concluded that the knowledge of the 
experimental vibrational distribution is important for a final conclusion. The present 
analysis should motivate further experimental studies, especially at variable photon 
energy, as it should become possible due to the current progress in the development of 
free-electron lasers. 
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